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A lens system for a periodic light-beam waveguide is proposed and 
analyzed in which gas is enclosed in a circular cylinder heated with a 
cos 2(f) temperature distribution. We show that this temperature distribution 
may be produced by cutting a cylindrical hole in the center of a square block 
which has two opposite sides of equal temperature above the ambient tem- 
perature, and two sides of a lower temperature. Heat conduction across 
the gas produces an index of refraction variation which, in two orthogonal 
azimuthal planes, increases or decreases as the radius squared. The effect 
of thermal convection is analyzed by solving the governing equations as an 
expansion in powers of the Rayleigh number; the solution 7 % eveals that 
convection effects can be made negligible over a practical range of lens 
parameters. The major attributes of the lens system are that only tempera- 
ture controls are required and the aberrations associated with thermal 
convection can be readily minimized. 

I. INTRODUCTION 

A gas lens system to transmit a light beam through a tube should 
have a favorable refractive index, negligible aberrations, and a simple 
construction. The favorable refractive index must be such that all 
light rays parallel to the tube axis, but of varying distances from that 
axis, converge at approximately the same point on the axis, the dis- 
tance being called the focal length. Within the paraxial ray approxi- 
mation it is easy to show that an r 2 variation of the refractive index 
has this property (see, for example, Refs. 1 and 2) . 

Berreman obtained a refractive index (which varied approximately 
as the square of the radius) by flowing a gas through a cold cylinder 
enclosing a warm helix aligned on the axis. 3 The interior of the 
helix has the desired refractive index. Marcuse and Miller simplified 
Berreman's lens by considering a cool gas flowing through a heated 
cylinder of uniform temperature (the Graetz problem). 1 - 2 
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In order to reduce the distortion resulting from spherical aberra- 
tions, Berreman built a counterflow arrangement composed of two 
back to back tubular lenses. 4 Marcuse calculated the principal sur- 
faces of a flow type lens noting that the one with the light beam 
parallel to the flow differs only slightly from that of the beam anti- 
parallel to the flow. 2 He then numerically calculated the fate of a 
beam as it passes through a large number of flow lenses and compared 
the results with those with a counterflow arrangement. 5 This arrange- 
ment decreased the distortion. Kaiser later found that this configura- 
tion also lessens the asymmetric distortion due to thermal convection. 6 

The major drawback to the flow-type lens is the need for control of 
the flow. Gu performed a compressible flow analysis and found that, 
as a result of the wall friction, choking could occur for the optimal 
flow rate in a few hundred meters. 7 This could be overcome only by 
further complexities in the system. 

A conduction-type lens was proposed by Suematsu, Iga, and Ito, 
in which they analyzed a configuration composed of hyperbolic, con- 
vex inward walls, two of which are at one temperature and the other 
opposing two at a lower temperature. 8 The concomitant temperature 
distribution varies as the square of the distance in the transverse di- 
rection. Then the refractive index bears the r 2 variation* in two 
orthogonal planes, being convergent in one and divergent in the other. 
This guadratic variation has two highly desirable characteristics. 
First, within the paraxial approximation, the focal length of every 
ray passing through a quadratic lens is independent of the radius, 
and hence the field reproduces itself after each period. 1 Second, 
Marcatili has shown that the eigenfunctions associated with a quad- 
ratic lens are Gaussian. Therefore a laser beam which is also Gaus- 
sian can be mode-matched to a waveguide consisting of quadratic 
lenses. This means that all the energy will remain in the launched 
mode; the only mode conversion that would take place is that result- 
ing from higher order variation, that is, aberrations. 

The advantage of the conduction lens is that only temperature con- 
trols are required since no gas flow is involved. However, thermal 
convection is present in this lens and although Suematsu, and others, 
observed a degradation of their lens at high temperature differences 
they did not analyze the thermal convection effects. 



* For negligible pressure changes the refractive index is virtually only a func- 
tion of the temperature; for small temperature variations the changes in refrac- 
tive index are directly proportional to and of opposite sign from the temperature 
changes. 
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We will consider a quadratic conduction-type lens, which is formed 
by imposing a cos 2</> temperature distribution on the wall of a cir- 
cular cylinder. For sufficiently small temperature variations the change 
in the refractive index is approximately quadratic and lensing action 
similar to that of Suematsu, and others, is obtained. The three cen- 
tral questions considered are: (?') What are the effects of thermal con- 
vection on the quadratic distribution? (u) How does one readily 
obtain a cos 2<f> wall temperature distribution? (Hi) What are the opti- 
cal properties of a waveguide consisting of these lenses? 

We show that the cos 2$ distribution can be achieved very simply by 
boring a circular hole in a square block in which two opposite sides bear 
a higher temperature than ambient and the other two bear a lower tem- 
perature. If sections of the above lens are placed in tandem, each con- 
secutive one rotated by 90 degrees, there then exists in one plane a 
scries of alternating, convergent-divergent lenses. In the perpendicular 
plane this series is, so to speak, 180 degrees out of phase. We may 
then use Miller's analysis of a sequence of alternating gradient 
lenses,* and determine criteria for the optical properties as a function 
of the parameters of the system. 

We study the effect of thermal convection by using a straightfor- 
ward perturbation analysis which is found to be in agreement with 
preliminary results of an experiment. We investigate the method of 
producing the wall temperature distribution by constructing an ap- 
proximate solution which reveals how the wall temperature distribu- 
tion can be established, as well as discuss the experimental program 
in progress and compare this lens and the other cited above. 

II. ANALYSIS 

2.1 Analysis of Thermal Convection 

Consider a circular cylinder with the geometry given in Fig. 1. The 
governing equations for the steady motion of the gas within the 
cylinder are: 

(i) continuity equation. 

V(pu)=0; (1) 

(ii) equation of motion, 

p(u-V)u + Vp - P g - mV'u = 0; (2) 

* Alternating gradient focusing in gas lens systems was first proposed by A. It. 
Hutson. 3 
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Fig. 1 — Geometry of the lens cylinder. 

{Hi) energy equation, 

pu>V(c p T) = kV' 2 T + u-Vp + /**■ 
Here, 



(3) 



p is the density, 

u the velocity, 

p the pressure, 

g the gravitational acceleration, 

n the viscosity, 

c p the specific heat at constant pressure, 

k the thermal conductivity, and 

<i> the dissipation function (associated with the frictional work). 

The boundary conditions at the cylinder surface are 



/ o 7 • (1 + |^ cos 20 



(4) 



and 



u(o, <*>) = 0, (5) 

where AT is the maximum excursion about the average wall tempera- 
ture, T . 

At this point we use the Boussinesq approximation which consists 
of two elements; the density changes are significant only in the body 
force term, and these changes are a function of temperature only. The 
latter element amounts to neglecting the product of the isothermal 
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compressibility, k, times the pressure change in comparison with the 
product of the volumetric expansivity, {3, times the temperature 
change. In other words, for p = p {p, T) 

= K(Ip - & dT, 
and the Boussinesq approximation requires the second term to be 
much larger than the first term but still small enough so that 

p = poll - {S(T - T.)], (6) 

where the subscript denotes conditions at the center of the cylinder 
in the absence of fluid motion. 

We nondimensionalize the variables in the hope that a perturbation 
scheme for a solution to our problem may be suggested. We define 

U-TTtS. x = X/a, B = ^^- (7) 

k/( Po c p a) Ai 

Since the density changes are considered important only in the body 
force term, equation (1) yields the incompressible continuity equa- 
tion, 

V-U = 0. (8) 

The pressure term can be eliminated from the equation of motion by 
taking the curl of equation (2) . The result of this operation leads us 
to define the velocity components in terms of the stream function \f/ so 
that in cylindrical coordinates we have 

fc-JSf. *--£ (9) 

The continuity equation, (8), is identically satisfied, and the equa- 
tion of motion becomes 

rV v + 'A d± --- t rV* = *(« * f - 8 -^ £)#, (io) 

* v V x r \dr d<f> d<j> drl \ dr r bty 

where 

a = 2£z the Prandtl number, and 
k 

= &Tgc vP 'tf the Rayleigh num ber. 
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Fig. 2 — Contour plot of the first approximation for the stream function with 
values of ^- (1) on indicated contours: 



A = 0.00001 

C = 0.00005 

D = 0.0001 

E = 0.0002 



F = 0.0005 
G = 0.001 
H = 0.002 



M = -0.002 
N = -0.001 
O = -0.0005 
P = -0.0002 



Q = -0.0001 
R = -0.00005 
T = -0.00001 



If the velocities are sufficiently small then the viscous dissipation 
can be neglected and the energy equation, (3), in terms of the new 
variable becomes 



r \dr d(f> d<t> arl 
The boundary conditions, (4) and (5) , become 

0(1, 0) = cos 20 
and 



(U) 



(12) 



2Ja.*-gfa.<)--o. 



(13) 



In the case of a small Rayleigh number it is fruitful to seek a solu- 
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tion in powers of A; 

* - X*" 1 + XV (i) + ••• (14) 

and 

= e w + \$ w + \*o m + ••• . (15) 

This expansion is valid in the limit, A -* 0, and an upper bound of A 
for the validity of the expansion will be obtained subsequently. 

When we insert equations (14) and (15) into (10) and (11), the 
coefficients of like powers of A must individually be set equal to zero 
for the equations to hold as A is varied. Beginning with the lowest 
order we obtain from equation (11) 

VV W = 0. (16) 

The solution to the equation, with the boundary condition given by 
equation (12) , is 

fl (O, (/-,0) = r cos 20. (17) 

Next, from equation (10) the lowest order contribution to the 
stream function is obtained from 

dr r d(j) 

with the boundary conditions given by equation (13). Inserting equa- 
tion (17) into equation (18) and expressing the biharmonic operator 
in cylindrical coordinates yield 

/of 4 , 2 d" 1 £_ , l_± 

\dr A + r d? ? dr 2 + r 3 dr 

-*7^r + ?^ + ?lz+? B)^ -* C ° S *■ (19) 
r d$ or r o<f> or r o<p r oq> 1 

The solution to this inhomogeneous biharmonic equation is 

<A ( >, <fi) = ^ [r 4 - °-r + l]r cos*. (20) 

Figure 2 is a contour plot of the stream function, equation (20) . 

Finally we wish to determine the perturbation on (o) . This will 
indicate the effect of thermal convection in distorting the lens and 
afford an estimation of the upper bound of the Rayleigh number. 
Again, from equation (11) we get 
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Fig. 3 — Contour plot of the (a) zeroth approximation for the temperature 
distribution with values of (o) on indicated contours: 



A = 0.1 = -R 
B = 0.2 = -S 
C = 0.3 = -T 



D = 0.4 = -TJ 
E = 0.5 = -V 
F = 0.6 = -W 



G = 0.7 = -X 
H = 0.8 = -Y 
I = 0.9 = -Z 



(b) first perturbation for the temperature distribution with values of (1 > on indicated 
contours : 



A = 0.0001 = -S 
B = 0.001 = -T 



C = 0.005 = -U 
D = 0.01 = -V 



E = 0.015 = -W 
F = 0.02 = -X 



G = 0.025 = -Y 
H = 0.028 = -Z 



(c) first approximation for the temperature distribution (X = 10 3 ) with values of 
0<°> -f X0 (1) on indicated contours: 



A = 0.1 = -R 
B = 0.2 = -S 
C = 0.3 = -T 



D = 0.4 = -U 
E = 0.5 = -V 
F = 0.6 = -W 



(i = 0.7 = -X 
H = 0.8 = -Y 
I = 0.9 = -Z 



„,„<„ l far ae™ 

V e ~ r\ d<f> Br 



(i) aa<°> 



d^" dd 



dr d<f> 



with the boundary condition 



U> (!,<*>) =0. 



Inserting equation (17) and (10) into equation (21) yields 



(21) 
(22) 
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+ 2(5r 4 - 6r 2 + l)r cos sin 20. (23) 
The solution which satisfies equation (22) is 

* (n (r,0) = [/.(l) - / 2 (D]|sin0 - [/,(1) + / 2 (l)]|sin30 

+ f,(r) sin cos 20 + / 2 (r) cos sin 20, (24) 



where 



and 



fM _ r* /32 48 ,\ 

hK) 96A379 2041 r I 



i (A - -!— (%■* 48 ( 31 ) J _i * -1_ 32 ( n ) r i o 
/a(l) = 4(96) [ 2l ~ -20ir r " 4r + ^79" ? + 2 

A numerical calculation reveals that the maximum value of (1) is 
approximately 3 X 10 -4 . Since (o) is bounded by unity, the expansion 
should be valid for Rayleigh numbers less than the order of 10 4 . Figures 
3a, b, and c show contour plots of (o) , (1) , and (o) + A0 (1) , respec- 
tively. In Fig. 3c, A = 10 3 to demonstrate the distortion possible. 

Experiments are being conducted to verify the foregoing results and 
to better understand thermal convection in other circumstances. Fig- 
ure 4 is a photograph of the streamlines made visible by the introduc- 
tion of cigarette smoke into a circular cylinder having a cos 20 
temperature distribution. The Rayleigh number is 575. Notice the re- 
semblance between this pattern and the contour plot of the preceding 
analytical results (Fig. 2). The slight shift upward of the smoke 
streamlines can be attributed to higher order terms in 6 and \f/. The 
steadiness of the observed flow supports our seeking time-independent 
solutions of the equations of motion. 

2.2 Establishing the cos 20 Wall Temperature Distribution 

If one imposes a linear temperature distribution across a slab by 
heating one face and cooling the other, and then if one drills a 
cylindrical hole parallel to the faces of the slab, it is well known that 
a temperature distribution varying as cos will appear on the wall 
of the cylindrical hole. Extending this to a square region with one pair 
of opposite faces heated and the other pair cooled one might presume 
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that a cos 2<f> temperature distribution would appear on a cylindrical 
hole cut in the center of the square. To determine the degree of ap- 
proximation of this presumption the heat conduction problem in a 
region bounded on the exterior by a square and on the interior by a 
circle is analyzed in the following paragraphs. Figure 5 shows the 
geometiy of the problem. 

The problem of a square with a hole in it cannot be solved exactly, 
as we show. An approximate solution could be sought in either cartesian 
or cylindrical coordinates. However, considering the problem in cy- 
lindrical coordinates allows one to compare the relative magnitude of 
the portion of the distribution, which varies with cos 2<f>, to that as- 
sociated with higher order terms. Secondly, the solution is more nearly 




Fig. 4 — Convective motion illuminated by cigarette smoke. 
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Fig. 5 — Geometry for conduction problem in Bolid cross section of gas lens. 

exact on the cylindrical hole if the approximate solution is sought in 
cylindrical coordinates. Furthermore, since the heating arrangement 
has a certain amount of symmetry, only a sector x/4 ^ <f> ^ x/2 need 
be considered. 

For steady two-dimensional conduction in a material having constant 
thermal conductivity the heat conduction equation becomes 



§11 _i_ ±§1 j_ 1.3-1 _ n 
-,2~r - ~r 2 ^ ,2 — "■ 
or r or r d<f> 



The boundary conditions 


are: 




at r sin <j> = b, 


6 = - 


at 


r = 1, 


dr 


at 


*-;■ 


= 0; 



at 



<t> = 



2 ' 



do 



= 0. 



(25) 

(26) 
(27) 

(28) 

(29) 



Notice that = (T — T )/ATa,s before, where T is the tempera- 
ture excursion desired on the cylindrical hole. Consequently, © = 
( Ty, — T ) /AT where T w is the wall temperature. The insulated condition, 
(27), assumes that the heat lost to the gas in the hole is negligibly small 
compared with the heat conduction in the solid. This is reasonable as 
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long as fc. oli d » fc gaB .* Condition (29) results from the symmetry about 
x/2. The radius r is normalized with respect to the cylinder radius 
as in the Section 2.1. 

Assume a separable solution of the form 

= Z2(r)$(0), (30) 

so that 

r 2 R" + rR' - aR = 0, (31) 

$" + a 2 $ = 0. (32) 

The solution of equation (31 ) and (32) is 

B = A\r" + — J(Csina0 + cosae£). (33) 

The insulated condition (27) is satisfied if B = 1. To satisfy both 
conditions (28) and (29) simultaneously, C = and 

a = 2n, n = 1,3,5, ••• . (34) 

Therefore, 

6 n = A n (r 2n + iVcos 2nct>), n = 1, 3, 5, ■ • ■ , (35) 

where -A„ is determined to satisfy equation (26), that is, 

-0= V A n (^4^ + sin 2n <*>) cos2?i(^. (36) 

Because of equation (26) our problem in r is not a Sturm-Liouville 
system and we have no assurance that equation (36) will converge 
even if the A n 's could be determined in general. In what follows we 
determine the first few 4 n 's so that equation (36) is satisfied in 
two different senses as accurately as our needs dictate — collocation 
and minimization of the error in a least-squares sense. 10 

In the collocation method the error is made to vanish at, say, three 
particular points on the boundary r sin = 6. This gives us three 
simultaneous equations through which A t , A 3 , and A 5 can be deter- 
mined. For two different sets of collocation points, the corresponding 
coefficients are listed in Table I for the ratio of the side length to the 



* The k for most plastics is a factor of 10 greater than that for air. For formed 
plastics taoiid «* fcgas and the behavior of at r =1 can be assessed from the solution 
for conduction in a square with two sides at and two sides at — 0. 
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Table I — Coefficients for Approximate Solution by Collocation 





Collocation points 


= 60°, 75°, and 90 






6=2 


6=4 


6=6 


ai 

a 3 


1.03531 

-0.11082 

0.01083 


1.08889 

-0.10469 

0.01155 


1.09191 

-0.10434 
0.01159 




Collocation points 


= 50°, 70 °, and 90 






6=2 


6=4 


6=6 


ai 
a 6 


1.04521 

-0.14096 

0.03046 


1.09927 

-0.13457 

0.03101 


1.10232 

-0.13421 

0.03104 



cylinder diameter, b = 2, 4, 6. These coefficients are normalized with 
respect to 0. Furthermore, some of the dependence on 6 is suppressed 
when the coefficients are denned as: 



a n = 



A n b 2n 




so that 



- „=tt 



vfe (r 2n + -2^j cos2n4>, 



(37) 



(38) 



Figure 6 contains a plot of 6{r = 6/sin <£, x/b) using both sets of col- 
location points. This illustrates the degree of approximation entailed 
at the outer boundary where —(0/©) should equal unity over ^ 
X/b < 1. 

The least-squares method requires that the mean square error over 
the boundary r = 6/sin <£, tt/4 ^ # ^ tt/2, be as small as possible. 
Defining 



= e - (-0); 



we then wish to minimize 



/: 



e d<$>. 



(39) 



(40) 



For convenience we take only the first two terms of equation (35) for 
6 and note that r 2n ^> l/r 2n close to the outer boundary so long as 
6^2. Consequently, 



8 = A x r 2 cos 2<f> + A 3 cos 60, 



(41a) 
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and 



t ^ Aib >^+ A3b *™<*> 



(41b) 



(-) 



1.3118 
b 2 



sin" " "~ sm <t> 
Inserting equation (41b) into integral (40), performing the integra- 
tion, and setting the derivative with respect to Ai and A 3 equal to 
zero we find that: 

( 2 , l\ , 0.1S05 /„ . l\ e , fjlo . 

[r + -2 ) cos 20 - , 6 lr 8 + ^J cos 60. (42) 

Figure 6 also has a plot of equation (42) evaluated at r = b/sin 0. 
Apparently the collocation method yields a much closer approxima- 
tion for heat conduction problems. (The square of the temperature 
has no particular physical meaning.) 

Returning to the collocation solution (Table I) the temperature 
distribution on the cylindrical wall, for b = 6, is given as 

0(1,0; 6) 2a, cos 20 2q 3 cos 60 2a 5 cos 100 
~@~ = 6 2 + 6 fl + 6'° 




Fig. 6 — Comparison of approximation at outer boundary in conduction problem. 
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& 0.06066 cos 20 - 0.0447 X 10" 4 cos 60 (43) 

+ 0.0383 X 10" 8 cos 100, 

and we see that the distribution varies as cos 20 within one part in 
10,000. For b = 2 the deviation is somewhat greater, being 

^'fr 2 ) = 0.5176 cos 20 - 0.00346 cos 60 + 0.02115 X 10~ 3 cos 100. 

(44) 

Similar results can be obtained from equation (42). 

Recall that this solution is for fc Bolid >s> k eaa . When /c Bolid tt fc g »„ the 
deviation of 0(1) from cos 20 can be evaluated from the analytical 
solution for a solid square two sides at and two other sides at —0, v. 11 
In doing this we found that the deviations from cos 20 are of the same 
order as those cited above. 

The power necessary to operate the lens can be readily found from 
integrating along the radial line at = 7r/4. The heat flow rate Q 
through one sector is given by 

Q = AT Zj A ,,^[(^)'' ~^] (45) 

which is nearly independent of b. In terms of the collocation coef- 
ficients, where the 6'"' term has been neglected, 

j^7 = 2ai -8a 3 + 32a 5 . (46) 

For b = 4 

Q = 3.3S fc6 AT = 3.38 k(T w - T B ). (47a) 

Since 

001,1-/2; 4) = -0.1360 = -0.136( - " ~ T " ) = 1, (47b) 

then for a AT of 1°C excursion T wM - T = 7.4° C. For a gas lens 
whose solid portion is made of polystyrene (fc = 0.1 W/m°C) the heat 
flow rate would be Q = 2.5 W/m for each sector; then the power re- 
quirement would be 10 W/m for 1°C AT across the lens. If a foamed 
polystyrene could be used the power requirement would be 3.5 W/m°C. 
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2.3 Optical Properties of the Conduction Lens 

With the effect of thermal convection in mind we now consider how 
to evaluate the optical properties of a gas lens characterized by the 
lowest order temperature distribution (that is, r 2 cos 2<f>). We wish 
to determine these properties as functions of A7 7 , cylinder radius a, 
and lens section length L; we are constrained by the requirement of 
minimizing A7 7 so that the distortion depicted in Fig. 3c shall be toler- 
able. In the following paragraphs we only write clown the relevant 
equation; we do not establish explicit design criteria. 

The system of lenses consists of a sequence of sections with each 
succeeding one rotated 90 degrees. Therefore, for any angle, <£ (sec 
Fig. 1), as one marches axially, the sections act alternately as diver- 
gent and convergent lenses. Since the temperature varies angularly, 
as well as radially, so does the refractive index; hence, in addition to 
the ray bending toward or away from the axis it will, in general, be 
twisted. However, at <£ = and tt/2 the refractive gradient has no 
angular gradient and, hence, rays originally in either of those planes 
remain there; they undergo convergent and divergent displacements 
alternately. All other rays have radial displacements intermediate 
to those at </> = 0, tt/2. 

The trajectories of the rays in the </> = 0, tt/2 planes may be calcu- 
lated analytically and turn out to be sinusoidal and exponential in the 
convergent and divergent sections, respectively. 

Although a numerical solution must be used for the other trajec- 
tories, some qualitative observations may be made. In the neighbor- 
hood of <}> = the angular component of the refractive index causes 
rays to be twisted away from that attitude, while near <£ = tt/2 rays 
are restored to that angular position. Therefore, as one moves down 
a section the density of rays tends to increase near <j> = tr/2 and to 
decrease near <f> = 0. 

In order to obtain the intensity of the beam through a lens section, 
the Helmholtz-type equation with the appropriate refractive index 
must be solved. This was done by Marcatili 12 for an asymmetrical 
but convergent-type refractive index." He established conditions for 
the stability of a lens system and calculated the focal length. 

For our present purposes there is no need for a detailed solution of 
the field equations but rather for the ray displacement, stability 
criterion, and focal length. Toward this end Miller's 9 analysis of the 



* Marcatili informed us that there is no basic reason why his analysis could not 
be extended to include divergent sections. 
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ray equation is applicable.* He obtained these quantities by solving 
the difference equations which govern the passage of the rays through 
the sequence of lenses. If we only consider the <f> = 0, tt/2 planes, then 
the sections act as alternating convergent and divergent lenses, with 
the rays remaining in their original planes; we may then apply Mil- 
ler's results. 

Miller obtained the ray displacement after the nth convergent and 
with divergent lens for an initially convergent and an initially diver- 
gent sequence. He also found the stability condition which keeps the 
ray trajectory bounded after an infinite number of lenses. This con- 
dition is 

< j < 2. (48) 

Here, to serve as an example, we only display the expression for the 
ray displacement after the nth convergent lens for an initially con- 
vergent lens: 

r„ = r a k\ cos (n 8 — <f>i) + r' Lk 2 sin n 8 (49) 

where r„ and r'„ are the initial displacement and slope, respectively, and 

—rb-wr ■ 2 



*, = 



i-r 



<f>i = I cos fc, -1 I, and fc 2 = 



1^ 

sin 8 



2/ 

(50) 



Furthermore, we must stipulate that the ray does not intersect the 
cylinder wall, that is, 



T i < l - < 5 » 



The relationship between the focal length and the refractive index 
may be obtained from Marcuse and Miller. 1 For a thin lens the focal 
length is given by* 



t We are indebted to Marcatili for several clarifying remarks on this subject. 

* A thin lens is one in which the principal surface generated by rays incident 
from the left coincides with that surface constructed by rays incident from tho 
right. Since there is no preferred direction with the conduction-type lens it is 
thin; the flow-type lenses cited in Section I may be approximately thin. 
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1~*k- (52) 

where p = 2ir/k, a is the wave length of the light, and &<f> is the dif- 
ference of the phase of a ray incident a distance r from and parallel 
to the axis after traveling a distance L, compared with a ray on the 
axis traveling the same distance. 

To calculate A<£ in terms of the refractive index, we invoke the 
paraxial approximation in which the rays are regarded as approxi- 
mately parallel to the axis. Then the required phases are easy to 
calculate, that is, 

0(r, z) ^ /3„ [ n(r) dx = o n(r)L (53) 

•'n 

and 

0(0, z) = 0„n(O)L. (54) 

The refractive index at <f> = is 



(55) 



T n„ — 1 

n(r, 0) = 1 + (n. - 1) ^r^ = 1 + ^2 

X + T n a 2 
,. AT r 2 

where n is the refractive index at the axis at temperature T„. 
Hence, 

AT r 2 

A0 = /3 o (n o -l)|r-^i W 

and the focal length is obtained from equation (52) : 

/ = I flS (57) 

7 2 (n. - 1) ATL 

independent of r. 

With the aid of equations (48), (49), (51), and (57) we may deter- 
mine the focal length and ray displacement as a function of the lens 
section and radius and temperature excursion. For a complete dis- 
cussion of the foregoing subject, see Ref . 9. 

To establish precise design criteria the foregoing equations must be 
solved on a computer. However, for illustrative purposes and as one 
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aspect of the problem we shall make use of some of Miller's simplified 
expressions valid in certain limits. 9 

If we use equation (49) with the value of the section length to 
focal length ratio which yields the smallest value of the maximum ray 
displacement and, furthermore, insure that the rays do not intersect 
the wall, AT obtained is unacceptable for three major reasons (i) the 
power requirement is excessive, (it) the moderate Rayleigh number 
will cause appreciable distortion, and {Hi) the temperature excursion 
is sufficiently large so that section end effects may be significant. 

In order to overcome these objections we now examine the case of 
weak focusing, that is, 2//L ^> 1. We consider the initial conditions 
such that 

r„ « r'J. (58) 

(The opposite inequality for weak focusing yields a trivial design 
problem since it does not involve the focal length.) From Miller the 
maximum ray radius, r mnx , is 9 

r m « = 2K ■ (59) 

To insure that the ray does not intersect the wall we have 

| i *£■ (60) 

Inserting equations (57) and (59) into equation (60) yields 

a T 

1 > - ±-£ -./ (Ql) 

= L(n - 1) AT ° ' K } 

As an example we use the following values, where air is the medium 
of the lensing action 

T = 290 °K, 

n„ - 1 = 0.295 X 10" 3 , 

and 

X (Rayleigh number) = 9.15 X 10 7 A TV 

with AT in degrees Celsius and a in meters. In addition, let r' Q = 10~ 4 , 
a = 3 X 10~ J m, and L = 0.5 m. Then from inequality (61) we obtain 

AT ^ 0.59°C. (62) 

Hence, X = 2.9. Consequently, 
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| \B U) | < 10 '•' « | d m U . 

It should be borne in mind that with the above values of a and L, the 
weak focusing limit is satisfied for AT < 5°C. In addition, to satisfy 
inequality (58) with the foregoing values we must have 

r„ « r'J = 1.5 X 10" 3 m (63) 

which is easy to satisfy. 

Using equation (47a), we obtain for the heat flow rate through one 
sector 

(J - 1.5 W/m 

which results in a power requirement of 6.0 W/m. The required ex- 
terior wall temperature is calculated from equation (47b) as T w = 
T + 4.4°C. Therefore, in the limit of weak focusing the temperature 
excursion is sufficiently small to make the lens system promising. 

Considering the flow-type lens of Marcuse and Miller to have the 
same characteristics as the above conduction type lens, we calculate 
the power expended at optimum flow rate to be 1.14 W. 1 Hence, the 
lens proposed here requires somewhat more power for heating than 
those previously investigated. However, the flow-type lens also re- 
quires power to drive the gas. 

Since the input beam will be more complicated then was assumed 
above, the foregoing calculation is very cursory. However, the rea- 
sonable magnitudes of a and L together with the small Rayleigh num- 
ber lend encouragement to a more detailed analysis. 

III. CONCLUSIONS AND RECOMMENDATIONS 

The conduction-type lens proposed here is found to be feasible on 
the basis of negligible distortion resulting from thermal convection 
and reasonable power requirements to maintain the desired tempera- 
ture distribution. Although the lens design illustrated was predicated 
on the weak focusing limit a wider range of parameters can be found 
by using Miller's complete expression. 9 

The effect of thermal convection was calculated from a two di- 
mensional analysis, which is certainly valid away from the ends of 
the section since a/L <& 1. For the temperature excursion required 
and the lens illustrated, the convection effect was found to be negligi- 
ble. However, at the interface between the sections, the axial tempera- 
ture gradients could be large depending on the spacing left between 
sections. Axial gradients were present in the experiments of Suematsu 
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and others for their hyperbolic shaped conduction-type lens system. 8 
They found that no significant aberration existed as long as AT < 
42/a 1 - 34 (a in millimeters) so that the effects of the axial gradients 
must have been insignificant. 

The analyses presented indicate that a system of conduction-type 
lenses might be practical for an alternating gradient light-beam wave- 
guide. Such a system would require straight square rods with a cylin- 
drical hole. Two sides of the rod would be heated while the other two 
would be held at a uniform and constant temperature. This could be 
done by attaching aluminum fins which project into a constant tem- 
perature heat sink to the cooled sides. Such a heat sink is available 
for buried systems since, at depths greater than about five feet, the 
surface temperature changes are virtually damped out. Therefore, 
cooing is not required. 

The hole in the rod would be of the order of 6 mm in diameter and 
the exterior could be as small as 2.4 cm across a face. Larger hole 
dimensions could be used but, for the same size beam and lensing ac- 
tion, the temperature difference and power requirement would in- 
crease proportionately. 

After only a preliminary design analysis, where the simplest of 
Miller's expressions have been used, parameters have been obtained 
in the weak focusing limit which yield a power consumption some- 
what greater than but of the same order of magnitude as flow-type 
gas lenses. 9 Additional investigations are, of course, necessary. The 
distortion of a gaussian beam as it is launched through a lens system 
should be numerically calculated (similar to Marcuse's study for the 
flow-type lens. 5 ) The effect of the axial gradients that will be present 
at the interface between two lens sections will have to be assessed 
through experimental measurements of the optical performance of 
such a lens system. 
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